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Abstract. 

We have performed narrow-band imaging observations with the Very Large Telescope, aimed at detecting the 
Lyman continuum (LC) flux escaping from galaxies at z ~ 3. We do not find any significant LC flux from our 
sample of two galaxies in the Hubble Deep Field South, at z = 3.170 and 3.275. The corresponding lower limits 
on the F1400/F900 flux density (per Hz) ratio are 15.6 and 10.2 (3-er confidence level). After correction for the 
intergalactic hydrogen absorption, the resulting limits on the relative escape fraction of the LC are compared 
with those obtained by different approaches, at similar or lower redshifts. One of our two objects has a relative 
escape fraction lower than the detection reported by Steidel et al. in a composite spectrum of z ~ 3 galaxies. 
A larger number of objects is required to reach a significant conclusion. Our comparison shows the potential of 
narrow-band imaging for obtaining the best limit on the relative escape fraction at z ~ 3. Stacking a significant 
number of galaxies observed through a narrow-band filter would provide constraint on the galactic contribution 
to the cosmic reionization. 

Key words, cosmology: observations — diffuse radiation — intergalactic medium — galaxies: photometry - 
ultraviolet: galaxies 

1. Introduction would gain additional support from direct measurements 

of the amount of Lyman continuum (LC) radiation re- 
The observations of the finishing and the beginning i ease d by galaxies into the intergalactic medium (IGM). 
reionization o f intergalactic hydrogen a t z ~ 6 and 20 

(Bec ker et ail l200ltlK~ogut et all 120031 co mbined with T _. _ il T ._ L . , 

T. T~" . ... . . , , ~~. 7v \ r / . ~~ 71 Direct measurements of the LC radiation from galax- 

the later reionization ol intergalactic He 11 (Kriss e t all . , . „ . T , . , ° 

|,.„„ ,1 rr, — — 71 \ nnn .1 „ i . les are, however, very challenging. In the local universe, 

2001; Zheng et al.. 2004) at z ~ 3 are now placing con- . . . . . T „ r 

, , , , ,. . , , they require space instruments sensitive at short U V wave- 

stramts on the nature and the evolution ol the ultra-violet , „ ,. . r -11 

/TT , T -, , n 1 t j- a 1 j -i , • r lengths as well as galaxies slightly redshiltcd to avoid the 

(UV) background radiation. A dominant contribution of j i- u • ^ 



star-forming galaxies is suggested at z > 3 before the den- 
sity of luminous quasars become significant. This picture 



neutral hydrogen opacity (continuum and line absorption) 
in the Galaxy. Observations of nearby star-forming galax- 
ies with the Hopkins Ultraviolet Telescope (HUT) and the 

~ ; ~ ~~ ; ~~ , -,, T Far Ultraviolet Spectroscopic Explorer (FUSE) have led so 

bend offprint requests to: A. K. moue , tp.pt/-i/ ■ 

* t, , , ,. „ . , . — c far only to upper limits on the fraction of LC (or emitted 

Based on observations collected at the European Southern „ J v 

Observatory, Paranal (Chile); Proposal No.: 71.A-0418(A,B). 900 A ) Photons that escape; the so-called escape fraction 
** Postdoctoral Fellow of the JSPS for Research Abroad. ( f mr ) is typically less than about 10% llLeitherer et all 



Present address: Okamaya Astrophysical Observatory, ll995tlHurwitz et allll997tlDeharveng et alll200l|h This is 



Honsho, Kamogata-cho, Asakuchi-gun, Okayama 719-0232, consistent with the values derived from the near-blackness 
Japan of interstellar lines tracing the H I gas as observed with 
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FUSE in five other starburst galaxies IjHeckman et all 
1200 ll) . 

At high redshifts, the increasing IGM opacity (absorp- 
tion of the ionizing radiation by the neutral hydrogen), not 
to speak of galaxies becoming fainter, makes the observa- 
tions difficult. Nevertheless, large gro und-based telescopes 
can enter the competition at z <; 3. ISteidel et al.l l)200l[) 
(hereafter SOI) derived a 1500A/900A observed flux den- 
sity 1 ratio, (i ? i5oo/^90o)obs = 17.7 ±3.8 from a composite 
spectrum of 29 Lyman break galaxies (LBGs) at a mean 
redshift of 3.4 with Keck/LRIS. By comparison with mod- 
els of the UV spectral energy distribution of star-forming 
galaxies, the flux density ratio, corrected for the IGM 
opacity, leads to the fraction of escaping LC (900A) pho- 
tons relative to the fraction of escaping non-ionizing UV 
(1500A) photons. This is called the relative escape frac- 
tion (/e S c,rel) and is different from the first definition of 
the escape fraction (/ esc ) used above for nearby galaxies. 
SOI interpreted their results as implying / eS c,rel <<S 50%. 

However, not all z ~ 3 galaxies emit a significant LC 
radiation. Using the FORS2 spectrograph on the VLT, 
iGiallongo et alJ |2002) (hereafter G02) obtained \-a lower 
limits on (-Fi5oo/^90o)obs four times larger than t he value 
of SOI for two bright LBGs. lHeckman et all l)200l|) also de- 
duced a very low escape fraction for a gravitational lenscd 
galaxies, MS 1512-cB58 at z = 2.7 from the detailed anal- 
ysis of interstellar absorption lines. 

Observations have not been only spectroscopic. 
Imaging of galaxies at z ~ 1 with the FUV solar-blind 
detector of Space Telescope Imaging Spectrograph (STIS) 
have also pro vided con s traint s on the flux below the 
Lyman li mit, llFergusonl | 200l| iMalkan et ail |2003). In 
particular. iMalkan et all l|2003() (hereafter M03) have ob- 
tained lower limits of (-Fi5oo/-FVoo)obs ^ 150-1000 (1-er), 
implying much lower LC escape fraction than in the galax- 
ies of S01. iFernandez-Soto et all l)2003|) (hereafter FS03) 
have used the deep U^oo images of the Hubble Deep Field 
North (HDFN) and reported an average LC escape frac- 
tion of no more than 4% for 27 galaxies at redshifts 
1.9 < z < 3.5. As a number of galaxies in their sample 
have their t/300 fluxes contaminated by nonionizing UV 
photons, their analysis is based on models and indirect. 

These conflicting results, as well as the difficulties of 
observations at high redshift, have led us to search for pos- 
sibilities of improvement. We have first examined how the 
detection advantage of imaging over spectroscopy is actu- 
ally working in the specific context of measuring faint LC 
radiation. On one hand, spectroscopy allows us to measure 
the LC flux close to the redshifted Lyman limit, where 
the average IGM opacity is not yet as large as it would be 
at shorter wavelength be cause of the Lyman valley (e.g., 
iMoller fc ,Iakobsenlll990h . On the other hand, broad-band 
measurements require a very large correction for the IGM 
opacity as we show in the current paper, because the ef- 
fective wavelength becomes very short for objects selected 

1 In this paper, all flux and luminosity densities are given 
per unit frequency interval. 



at a redshift appropriate for avoiding contamination by 
nonionizing UV photons. The narrow-band imaging ap- 
pears, therefore, as a natural compromise between the 
lower opacity offered by spectroscopy and the depth of 
detection offered by broad-band imaging. 

We have, for the first time, attempted imaging observa- 
tions from the ground, using the VLT, aimed at measuring 
the LC radiation from galaxies at high redshifts through 
a narrow-band filter. The goal of this paper is twofold: 
(i) reporting the data and their interpretation in terms of 
the relative escape fraction, and (ii) validating the narrow- 
band photometric approach by quantitative comparisons. 

In Sect. 2, we describe the details of our strategy, the 
choice of the filter, the target selection based on photomet- 
ric redshifts, the spectroscopic observations for confirming 
redshifts and the imaging observations. Data reduction 
and results are described in Sect. 3. Their interpretation 
in terms of relative escape fraction are presented in Sect. 4. 
A comparison of the constraints reached with those from 
other observations is given in Sect. 5, where we also discuss 
the absolute escape fraction. The final section is devoted 
to our conclusions. 

We adopt a flat A-dominated cosmology with fl\ = 0.7 
and fl m — 0.3, and the current Hubble constant Hq = 70 
km s _1 Mpc -1 to estimate luminosities. 

2. Observations 

2.1. Filter selection 

A direct and unambiguous detection (or upper limit) of 
the LC escaping from galaxies requires a lower redshift 
limit of the target galaxies to push out all emission red- 
ward of the Lyman limit beyond the long wavelength cut- 
off of the filter transmission. The cut-off wavelength is 
defined practically at 1% transmission relative to the peak. 
Since the IGM opacity rapidly increases with the source 
redshift (e.g., see Fig. A. 2), the limiting redshift should be 
set as low as possible. The Oil + 44 narrow-band filter, 
whose long wavelength cut-off is the smallest among all 
filters of FORS on VLT, was selected for our observation. 
Accordingly, galaxies with z > 3.18 must be selected. We 
note here that the long wavelength cut-off of this filter is 
smaller than that of the U 300 of WFPC2 on the HST as 
shown in Fig. 1. Indeed, the redshift limit for the U300 
would be z > 3.33. 

2.2. Observed field and target galaxies 

As the redshift surveys in the southern celestial sphere 
had not produced a large number of galaxies with spec- 
troscopic redshift (z sp ) larger than 2, we had to select 
galaxies based on photo metric redshif t (z n h ) in the first 
phase of our observations. lLabbe et all l)2003|) made a cat- 
alog of photometric redsh ift of galaxies in the Hubble 
Deep Field South (HDFS; IWilliams et all l2000l) by us- 
ing their ultra-deep near-infrared images obtained with 
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Table 1. Details of the imaging observations. 



2400 2800 3200 3600 4000 
wavelength [A] 

Fig. 1. Transmission curves of the FORS Oil + 44 filter 
(solid) and the HST U 300 (dashed). 



VLT/ISAAC and HST/WFPC2 optical images. Thus, we 
selected the HDFS/WFPC2 field as our target field. 

There are 24 galaxies with Jgi4 < 25.0 mag(AB) 2 and 
Zph > 2.5 in the field. We set this lower limit of z p h, taking 
into account the 2-a uncertainty of z n h for z ~ 3 (a ~ 
0.08(1 + z) noted in lLabbe etaP l2003|l . Eight of them 
have z p h larger than 3.18. On the other hand, six out of 
24 have z sp , but all of them are less than 3.18. We note 
that all the known z 

Zph- 



sp 



are smaller than the corresponding 



An accurate redshift is essential to avoid contamina- 
tion from nonionizing photons. We need a redshift accu- 
racy of the order of 0.01. Since the accuracy of z p h is not 
sufficient, we performed spectroscopic observations with 
FORS2 on VLT using the MXU mode. A full description 
of the observation will be published elsewhere (Iwata et 
al. in preparation). Here, we summarize the results briefly. 

Because of the limitation of the slit-let configuration, 
we were able to take spectra of only 15 galaxies out of 
the above 24 galaxies. Five out of the 15 galaxies have 
known z sp and 7 of them have z p h > 3.18. The choice of 
galaxies was made with a priority given to more luminous 
galaxies and galaxies with z p h > 3.18. Additionally, we 
took spectra of galaxies with z p h > 2.5 but Isi4 > 25.0 
mag(AB), for which we could configure the slit-let. Among 
these galaxies, one has z sp previously measured. We ended 
up with obtaining 13 z sp : seven of them are new mea- 
surements and the remaining six are confirmations of the 
previous results. 

Among the new redshifts, we got only one galaxy 
with a high enough redshift, z sp = 3.275 (HDFS 1825), 
in contrast to the expectation from photometric red- 
shifts. Indeed, for the 13 galaxies with z sp , we always find 



Dates (year 2003) 
Instrument 
Collimator 
CCD 

Pixel scale 
Field of view 
Filter 

Central wavelength 
Transmission FWHM 
Seeing 

Effective exposure time 



June 29, 30; July 5, 28; August 1, 28 
FORS1 on VLT 
standard resolution 
2048 x 2046 pixels 
0."2/pixel 
6.'8 x 6/8 
Oil + 44 
371.7 nm 
7.3 nm 
~ 0."8 
40,636 sec 



z S p < %>h, again. There is another galaxy with z sp = 3.170 
(HDFS 85) which was confirmed by our spectroscopy. This 
redshift is slightly smaller than the limiting redshift, how- 
ever, the transmission efficiency of the Oil + 44 filter at 
the Lyman limit of the galaxy is as small as a few per- 
cent of the peak efficiency (see Fig. 1). The contamination 
of nonionizing photons is small enough, especially in the 
context of upper limit measurements. Therefore, our final 
sample consists of the two galaxies. We have no other ex- 
planations for this small number of target galaxies than a 
systematic error in the estimation of z p h- 

2.3. Imaging observations 

The HDFS field was observed through the Oil + 44 
narrow-band filter with FORS1 and the TK2048EB4- 
1 detector chip in service mode from 29 June 2003 to 
28 August 2003. The standard resolution collimator was 
adopted, giving a pixel scale of 0."2 per pixel. The field 
of view is 6/8 x 6/8 and covers the whole area of the 
HDFS/WFPC2 field. Each image was taken with a small 
dithering (typically 10") and a typical exposure time of 
1,080 sec. The seeing size of each image was 0. "39-1. "23, 
typically ~ 0."8. We secured 38 images of the field. The 
effective exposure time is 40,636 sec (i.e. about llh). Table 
1 shows a summary of the observations. 

3. Data reduction and results 

3.1. Final image 

The image data reduction was carried out in a standard 
manner, using IRAF 3 . The bias subtraction was made us- 
ing over-scanned regions. Normalized twilight sky frames 
in each observing night were used for flat fielding. We se- 
lected 25 stars in the field and used them to register the 
38 frames. After registration, the rms of residual shifts 
was 0.03 pixel. These stars were also used for an airmass 
correction. We found a slope of magnitude dependence on 
airmass of —0.45, and we corrected observed counts as 



2 The is o-photal magnitude in the HDFS WFPC2 catalog 
version 2 dCasertano et all I2000T) 



3 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 




HDFS1825 




Fig. 2. Close-u p of sample ga l axies through the Oil + 44 filter (VLT/FORS, this work), and the isi4 filter 
CHST/WFPC2. IWilliams et al. L EoOO^ . The FOVs are 20" x 20". For the object 85, the image center is 6" shifted 
from the object position, because it is close to an edge of Igi4 image. 



it would have been observed at the zenith. The stability 
of the observing conditions is confirmed by measuring the 
corrected counts of these stars; the rms errors in counts are 
less than 0.05 mag for 19 stars brighter than 22 mag(AB) 
through the Oil + 44 filter among the above 25 stars. 
Then, the 1RAF task IMCOMBINE was used to sum the 
frames. We took averages of each pixel adopting a 3-er 
clipping. FWHM of stellar objects in the final image is 
~ 1". Fig. 2 shows the close-up images of the two sample 
galaxies through the Oil + 44 filter and through the I%u 
filter of HST/WFPC2. Neither galaxy seems to be seen 
through the Oil + 44 filter. 

3.2. Photometry 

Three standard stars, Feige 110, G 93-48, and LTT 9491, 
were also observed through the same filter. These frames 



were processed in the same way as the HDFS frames. The 
absolute AB magnitude of these stars were calculated from 
the spectral table provided on ESO web page. We found 
that the airmass dependence for the standard stars was 
—0.58 which is somewhat larger than that from stars in the 
HDFS frames. We derived the photometric zero point for 
the Oil + 44 filter at the zenith as 22.98 from the airmass 
slope of the standard stars. The zero point changes less 
than 0.02 mag if we adopt the airmass slope —0.45 which 
we used for correction of the HDFS frames. 

To reduce background fluctuations, we apply a 
Gaussian smoothing to the final image adapted to the ex- 
pected size of the sample galaxies through the Oil + 44 
filter. The size through the HST/WFPC2 B 450 filter is a 
reasonable approximation because the central wavelength 
is close to that of the Oil + 44 filter and we are seeing 
the light from massive stars in both filters. Diameters of 
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Table 2. Photometric properties of sample galaxies. 



ID a 



85 



1825 



a (J2000) 
8 (J2000) 
redshift 

Oil + 44 6 (nJy) 
U 300 b (nJy) 
-B4so c (nJy) 
Veoe c (nJy) 
J 8 i4 C (nJy) 



22 32 46.91 
-60 31 46.9 
3.170 

< 44.88 

< 101.8 
392.08 ± 12.02 
698.68 ± 8.86 
859.63 ± 17.82 



22 32 52.03 
-60 33 42.6 
3.275 

< 39.81 

< 42.09 
181.98 ±6.51 
404.91 ± 4.26 
504.49 ± 7.58 



a HDFS WFPC2 catalog version 2 dCasertano et al.ll2000lh 
b 3-<t upper limits. C/300 values of HDFS 85 and 1825 are 
reported in the HDFS WFPC2 catalog version 2. As these val- 
ues are similar to the 1-a uncertainty (HDFS 85) or negative 
(HDFS 1825), we have retained 3-er upper limits. Note that 
HDFS 85 is located near the edge of the WFPC2 image 
c With 1-a errors 



0."72 (HDFS 85) and 0."41 (HDFS 1825) are estimated 
from the half-light radii of the galaxies through the -B450 
filter reported in the H DFS WFPC2 catalog version 2 
ijCasertano et all l200fj|) . We adopt 2.5 x 2.5 pix as the 
smoothing scale. The FWHM of the PSF in the smoothed 
image is 8.33 pixel (= l."7). 

None of the sample galaxies was detected at a 3-cr level 
in the smoothed Oil + 44 image. The upper limit on their 
flux density has been based on the smoothed PSF size 
because they are smaller than the PSF. The measured 
rms fluctuations within a square box of l."7 x I. "7 in 
the smoothed image around the sample galaxies translate 
into 3-cr limiting magnitudes of 27.37 mag(AB) and 27.40 
mag(AB) for HDFS 85 and 1825, respectively. 

In Table 2, we summarize the photometric measure- 
ments of the two sample galaxies. Photometric data with 
H ST are taken from the HDFS WFPC2 catalog version 
2 l|Casertano et all l2000|) . In this table, all upper limits 
are 3-<r level. Our limiting magnitude is comparable with 
that of the HDFS in the U300 for the size of our sample 
galaxies, although the HDF S limiting magnitude fo r point 
source is deeper than ours i Casertano et alll200f)lh 



4. Lyman continuum escape fraction 

From the upper limits obtained above on the LC flux of 
the two galaxies at z = 3.170 and z — 3.275, we try here to 
estimate the escape fraction of LC photons. As mentioned 
in the introduction, the term of escape fraction has been 
used in at least two ways. We first clarify the definition of 
the escape fractions. 



4.1. Definition of escape fractions 

We define two escape fractions. One is 



j out 
T int 



exp(- 



TSMn 



(1) 



where is the intrinsic LC luminosity density (per Hz) 
of a galaxy, L£c 1& tne LC luminosity density just outside 
of the galaxy (not the observed one, see below), and t^q 1 
is the opacity of the interstellar medium (ISM) for LC 
photons in the galaxy. The other is 



csc.rcl 



T out 

T hit 

L uv 



— /esc exp(r uv ) , 



(2) 



where L^y, 



T OUt 

^TJV> 



and are the intrinsic luminosity 

density, outside luminosity density, and the ISM opacity 
of the nonionizing UV photons, respectively. Hereafter, 
the former is called the absolute escape fraction, and the 
latter is called the relative escape fraction. 

The IGM opacity should be taken into account for pho- 
tons with a wavelength shorter than the Lya line at the 
source rest-frame. Namely, the observed luminosity den- 
sity becomes L^ bs = i^ ut exp(— t{ gm ) for the rest-frame 
wavelength A < 1216 A, and i° bs = i° ut otherwise if the 
IGM dust extinction is negligible. Indeed, based on obser- 
vations of di stant supernovae and on the thermal history 
of the IGM. Ilnoue fc Kamaval l|2004l) have recently shown 
that the IGM dust extinction against z ~ 3 sources should 
be less than 1 mag in the observer's B-band. Taking a 
wavelength longer than the Lya line for the UV wave- 
length (e.g., ~ 1500 A), we have 



esc, rcl 



(Luy/iyLc)int 
obs 



exp(r^ M ). 



(3) 



where we have replaced the observed luminosity density 
ratio into the observed flux density ratio. Hence, the 
relative escape fraction can be estimated from the ob- 
served UV-to-LC flux density ratio if we know the intrinsic 
UV-to-LC luminosity density ratio and the IGM opacity. 
Moreover, the absolute escape fraction can be estimated 
from the relative escape fraction via equation (2) if we 
know the ISM opacity for nonionizing UV photons. 

4.2. Intrinsic luminosity density ratios 

The intrinsic UV-to-LC luminosity density ratio is still 
very uncertain observationally. As noted by S01, we must 
rely exclus ively on models. Here we adopt the Star burst 
99 model l(Leitherer et all Il999lh From this model, we 
obtain (Li5oo/£90o)int = 1-5-5.5 as shown in Fig. 3, as- 
suming a constant star formation rate, the Salpeter initial 
mass function with the mass range of 0.1-100 M Q , and 
the metallicity of 0.001-0.02 (the solar value is 0.02). The 
luminosity density ratios mainly depend on the duration 
since the onset of star formation; the ratio starts from a 
small value, monotonically increases with time, and satu- 
rates at a larger value after several 10 s yrs. In this paper, 
we adopt (£i5oo/£90o)int — 3.0 according to the value 
adopted by SOI 4 



4 S01 adopted the value as the Lyman discontinuity, and 
assumed the flatness of the intrinsic spectrum between 912 
and 1500 A. 
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Table 3. Adopted intrinsic ratios of luminosity densities 
relative to 1500 A. 



X (Ll50o/L x )i n t a 

IAJ 

700 4.0 

900 3.0 

1400 1.0 

1700 1.1 

2000 1.2 

2500 1.4 

a The unit of luminosity densities is erg s _1 Hz 




0l I Li LL_I i L 



500 1000 1500 2000 2500 3000 
rest-frame wavelength [A] 

Fig. 3. Intrinsic ratio of luminosity density at 1500 A to 
that at the wavelength along horizontal axis. Thin dotted 
curves show 16 cases of the combination of 4 metallicities 
(0.001, 0.004, 0.008, and 0.02) and 4 ages (0 yr, 10 Myr, 
100 Myr, and 1 Gyr), assuming a constant star formation 
and the Salpeter initial mass function with 0.1-100 M Q . 
Open squares are fiducial ratios adopted in the text at 7 
wavelengths marked by thick dotted vertical lines. Unit of 
luminosity densities is erg s _1 Hz -1 . 

Since the UV and LC wavelengths actually depend on 
the filter adopted and the source redshift, we have to ad- 
just the intrinsic UV-to-LC luminosity density ratios to 
the correct wavelengths. Some of the ratios, relative to 
1500 A, are displayed in Fig. 3 as open squares and sum- 
marized in Table 3. This table will be used to calculate the 
UV-to-LC luminosity density ratio appropriate for each 
individual galaxy in Sects. 4.4 and 5.1. As illustrated in 
Fig. 3, the uncertainties on the luminosity density ratios 
are 20-50%. 

4.3. IGM opacity 

An average opacity of the IGM can be estimated from 
the number distribution functions of the Lyman a forest 



and denser H I absorbing clouds fe.g.. lMadaulll995]) . The 
number of H I clouds per unit interval of neutral hydrogen 
column density JVhi and per unit redshift z i nterval can 
be exp ressed as d 2 N '/ 'dzdN\n oc (1 + zVN^f . iKim et all 
(2002) have shown that the observed Lyman a forest data 
can be fitted with (3 = 1.5, and 7 = 2.5 for z > 1.1 and 
7 = 0.2 for z < 1 .1, by combining their high-z data with 
the low- z data of lWevmann et all l)l998[) . 

In contrast to previous approaches fe.g.. lMadauL ll995'). 
we have shown in Appendix that not only the column 
density distribution but also the number evolution of the 
denser clouds (essentially, Lyman limit systems) can be 
reasonably described by the sa me function as t hat ob- 
tained for the Lyman a forest bv lKim et alJ (2002). Based 
on this distribution function, we calculate an average IGM 
opacity, (r^ GM (zs)), for a source with a redshift, zg, at 
an observed wavelength, A bs by an analytic approxima- 
tion formula presented in Appendix as equations (A. 4-6). 
As seen in Appendix, our opacity is quantitatively very 
consistent with the value measured and adopted by S01, 
an d shows a reasonable a greement with the opacity model 
of iBershadv et al] l)l999h . Finally, we integrate the IGM 
opacity over the filter transmission when we make a correc- 
tion for the IGM absorption in order to derive the escape 
fraction from the photometric data. In the integration, we 
always assume a fiat spectrum because we do not know 
the spectral shape of LC exactly. 

4.4. Estimation of relative escape fraction 

By combining the upper limit of the flux density through 
the Oil + 44 filter with the measurement of the flux 
density through the HST's Vao6, whose wavelength cor- 
responds to the rest-frame ~ 1400 A for our sample, we 
obtain a lower limit on the observed ratio (Fi4oo/Fgoo)obs- 
This ratio has to be corrected for the Galactic extinction 
before estimating the relativ e escape fraction. For HDFS 
field, E{B~V) = 0.028 mag llSchlegel. et al!ll998h . which 
corresponds to the Galactic extinction of Ay 606 — Aon = 
—0.050 mag if we adopt the standard extinction law of our 
Galaxy (e.g.. lWhitteiIl2003h . Therefore, we obtain the 3-ct 
lower limits of the corrected flux density ratio F1400 / -F900 
of 15.57 and 10.17 for HDFS 85 (z = 3.16) and HDFS 
1825 (z — 3.27), respectively. 

Assuming (Li4oo/£90o)int = 3.0 as found from Table 3 
and mean IGM opacities integrated over the filter trans- 
mission, we obtain 3-cr upper limits on the relative escape 
fraction as 72% and 216% for HDFS 85 and 1825, respec- 
tively. As found in equation (2), the relative escape frac- 
tion is not restricted to less than 100%, because it is the 
escape fraction relative to that of nonionizing UV photons. 
The obtained values are summarized in the top 2 lines of 
Table 4. 

Since uncertain assumptions on the ISM dust extinc- 
tion are inevitable for the calculation, the estimation of 
the absolute escape fraction is left for the following sec- 
tion (§5.2), where we will discuss these issues in detail. 
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Table 4. Summary of related observations and escape fractions. 



Redshift range/method/ 


redshift 


(jFjJY j i^LC )obs a 






Auv 5 -^lc b 


Z/tjv c 




ref. e 


/esc™ 


galaxy name 








T%) 




(10 29 erg s- 1 Hz" 1 ) 


(mag) 




(%) 


(1) 


(2) 


(3) 


(4) 


(5) 


((>) 


(7) 


(8) 


(9) 


(10) 


z ~ 3 

Narrow-band photometry (this work) 


















nur o oo 


^ 1 70 


> 15.57 


1.32 


< 72 


1 


1.59 


1.6 


1 


< 17 


riij r o lozo 


97^ 


> 10.17 


1.99 


< 216 


1 


0.970 


1.9 


1 


< 38 


Spectroscopy 




















composite (29 galaxies) 


q /in 


17.7 ±3.8 


1.50 


76 ± 16 


2 


1.60 


0.26 




60 ± 13 


Uar ZZo j -rllO UZ 


*3 Q1 Q 


> 24 


1.41 


< 51 


2 


3.15 


3.8 


2 


< 1.5 


onnnn oaq T\fi 


9 Qfil 

z.yoi 


> 23 


1.04 


< 37 


2 


4.85 


2.3 


2 


< 4.5 


Broadband photometry 




















PTVQQ- 

r Li i yy. yo < 


o.oO i 


> 2.34 


4.63 




3 


0.116 








PTVQQ- 89^ 
r L i yy. oZO 


o.ooy 


> 11.9 


4.64 




3 


0.617 








r _l i yy. oz4 


Q A QQ 


> 7.89 


4.94 




3 


0.638 








L? V A "\ Q /H r\QTl/H rl /""*L rif 

XDI UdLlUdllU. UllU LiUIIiobl V 






















1 1 074 

1 .11) 1 4: 


> 187 


0.53 


< 3.0 


5 


0.967 


1.3 


3 


< 0.91 




1 ISIS 

_1_ .1010 


> 308 


0.59 


< 1.7 


6 


2.56 






< 1.7 


r^PT5Q If! 1 SS7 


1 9Q7Q 

l.Zo f U 


> 119 


0.63 


< 4.5 


6 


0.969 






< 4.5 


r^PT5Q in noQQ 


1 9Q1 Q 

i.zyiy 


> 228 


0.68 


< 2.5 


6 


2.28 






< 2.5 


ppp q m 11 fiS 
iio lU.llOo 


1 1 ^Q9 


> 163 


0.57 


< 3.1 


6 


0.802 






< 3.1 


T HQQO 1 n 9CS 
IjUodz 1U. Zoo 


i i ns 


> 102 


0.53 


< 5.5 


5 


0.859 






< 5.5 


HDFiw4 1002 1353 


1.221 


> 67.8 


0.62 


< 10.0 


4 


0.513 


2.6 


4 


< 0.91 


CFRS 14.0547 


1.160 


> 348 


0.57 


< 1.8 


4 


1.55 


2.2 


4 


< 0.24 


CFRS 14.0154 


1.1583 


> 192 


0.57 


< 3.4 


4 


0.610 


2.9 


4 


< 0.24 


SSA 22-16 


1.36 


> 165 


0.74 


< 4.6 


4 


1.05 


0.20 


4 


< 3.8 


CFRS 22.1153 


1.3118 


> 125 


0.54 


< 5.0 


4 


1.00 


1.6 


4 


< 1.1 


Nearby galaxies 




















Spectroscopy 




















Mrk 54 


0.0448 


> 111 




< 2.7 


2 


0.83 


0.92 


5 


< 1.2 


Mrk 496 


0.0293 


> 7.5 




< 40 


2 


0.13 


2.4 


6 


< 4.4 


Mrk 66 


0.0218 


> 7.6 




< 39 


2 


0.064 


1.2 


6 


< 13 


Mrk 1267 


0.0193 


> 7.8 




< 38 


2 


0.068 


3.5 


6 


< 1.5 


IRAS 08339+6517 


0.0187 


> 12 




< 25 


2 


0.34 


1.6 


6 


< 5.7 



a Significance level of lower/upper limits is 3-<r, but not included uncertainties of model parameters. Uncertainty for Steidel's 
composite spectrum means 1-cr error of observations. 

6 Wavelengths of ultraviolet (UV) and Lyman continuum (LC); 1: (Auv,A L c) = (1400A, 900A), 2: (1500A, 900A), 3: 
(1400A,700A), 4: (1700A, 700A), 5: (2000A, 700A), 6: (2500A, 700A). 

c UV luminosity density corrected for the Galactic extinction. 

d In the estimation, j4i6oo-/3 relation by Meurer et al.(1 999) is assumed, and Ai6oo values are converted to those at the 
appropriate UV wavelength by Calzetti's attenuation law dCalzetti et all l2000l) . For Steidel's composite spectrum, we have 
used the mean E(B — V) value reported in their paper. 

e Reference of UV slope; 1: Veoe-hu color in HS T/WFPC2 catalog v ersi on 2 iCasertano et afl.l200(ft : 2: |Pettini et"all <ll998T) : 
3: B - V color in M03; kUn-G color in M03; 5: lMeurer etaH <1999h : 6: iLeitherer et all feOOaf) 



5. Discussions 

5.1. Comparison with other observations of Lyman 
continuum 

In the following section, we compare the limits obtained 
in Sect. 4.4 for our two galaxies with those already pre- 
sented in the literature (see Sect. 1). To do this, we have 
compiled all the UV-to-LC flux density ratios published 
in refereed journals. The data are summarized in Table 4, 
where they are divided into several categories depending 
on the source redshift and the observational technique. 



For the broad-band U300 data, we restrict ourselves to the 
galaxies with z sp > 3.33 in order to avoid any contami- 
nation of nonionizing UV photons. Thus, we took three 
galaxies from the list of FS03. 5 

We first compare the observed UV-to-LC flux density 
ratio (F\jy / Flc, column 3 in Table 4) and the estimated 



5 There is one more galaxy (HDFN:[FLY99] 688) satisfying 
the redshift criterion in the list of FS03, but it is detecte d 
by Chandra as a X-ray source iHornschemeier et all 1200 ll) . 
Hence, we remove it from our discussions because it may be an 
AGN. 
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relative escape fraction (/ CS c,rci, column 5), using the in- 
trinsic luminosity density ratios (Sect. 4.2) and the IGM 
opacity model (Sect. 4.3, column 5), which are tailored to 
the specific wavelength and the filter bandpass (column 6) 
for each individual galaxy. 

5.1.1. z ~ 3 galaxies 

The column 3 in the top of Table 4 shows that our two 
observed lower limits of F\jy / F^q are comparable to those 
already obtained for galaxies at z ~ 3. G02 have reached 
constraints about two times better than ours, because 
their two objects are more than two times intrinsically 
brighter. Our objects, as those of SOI, are comparable 
in luminosity to the of z ~ 3 LB Gs (see column 7; 
LJz = 3) ~ 1 x 10 2 9 erg s- 1 Hz" 1 : ISteidel et~Hl 
lAHelherger Stei dell . 120001 . 

In terms of relative escape fraction (column 5), we add 
one more object (HDFS 85) to the two galaxies of G02 that 
have been reported with a relative escape fraction lower 
than the detection of SOI. Because of the small sample 
size, these lower limits remain consistent with the rela- 
tive escape fraction of SOI. Although we have applied a 
mean IGM opacity for simplicity, a line of sight with larger 
than average opacity cannot be ruled out for an individual 
galaxy (not to speak of an unusual object as the escape is 
probably a random phenomenon) . 

The broad-band U^qo photometry in the HDFN 
reaches similar lower limits on i*uv /Flc (it would have 
reached slightly better values at the same luminosity level 
as ours). This, however, does not result in interesting con- 
straints on the relative escape fraction because of the IGM 
opacity. The IGM opacity through the [/300 filter, with an 
effective wavelength of 700 A in the source rest-frame, is 
much larger than those at wavelengths closer to the rest 
900 A (column 4). Since meaningful constraints on the rel- 
ative escape fraction are not reached, this sample will be 
removed from discussions in terms of the absolute escape 
fraction presented below. 

The top of Table 4 is also interesting for a comparison 
of methods of observations and a validation of future ap- 
proaches. As expected, the narrow-band photometry can 
provide better constraints than spectroscopy. The better 
constraints of G02 are actually due to the high luminos- 
ity of the objects. If our sample had a similar luminosity 
to the galaxies of G02, say, L = 4 x 10 29 erg s _1 Hz -1 , 
we would have obtained Fuv/^lc ^ 40 (3-cr level). The 
broad-band photometry is disqualified by the IGM opac- 
ity contribution (or requires the help of models that lead 
to less direct constraints). The narrow-band photometry 
can provide, even for a single LBG, a lower limit on 
-fuv / Flc comparable with that measured in the compos- 
ite spectrum of 29 galaxies by SOI. Stacking a significant 
number of galaxies observed through a narrow-band fil- 
ter will be able to go much deeper. We can also average 
the IGM opacity against any unusual line of sight in the 
stacking process. Such an approach would allow a signif- 



icant comparison with the result of SOI and reveal which 
fraction of high-z galaxies have high relative LC escape 
fraction. 

5.1.2. z ~ 1 galaxies 

Thanks to the high sensitivity of the HST/STIS, M03 have 
obtained very good lower limits on the flux density ratio 
of 11 galaxies at z ~ 1, -Fuv /-Flc > 70-350 at 3-cr con- 
fidence level (column 3 in the middle of Table 4). With 
the correction for the average IGM opacity described in 
section 4.3, these ratios are converted to upper limits on 
the relative escape fraction of 2-10% (column 5). The lu- 
minosities of this sample are similar to L* of z ~ 3 LBGs 
(column 7). 

The average IGM opacity (column 4) is much smaller 
than those at high redshift, but not completely negligible 
as assumed by M03, if we consider the Lyman limit sys- 
tems (LLSs) which dominate the IGM opacity for z ~ 1 
galaxies. The rarity of LLSs favors a statistical treatment 
rather than the average opacity adopted above. Based on 
the number distribution function of the IGM clouds as- 
sumed here (see Appendix), the expected number of LLSs 
within the wavelength range observed by M03 (1300 A 
^> A bB 1900 A) is about 0.3. Roughly speaking, for 
one-third of the objects of M03, the LLSs would loosen 
the constraints obtained under the assumption of no IGM 
opacity. In this sense, the number of galaxies in M03 is 
large enough to conclude a very small relative escape frac- 
tion for their objects. 

5.1.3. Nearby galaxies 

Although the upper limits on the LC from nearby galax- 
ies have been estimated from thei r Ha fluxes dire c tly in 
terms of absolute es c ape fraction llLeitherer et all Il995t 
iHurwitz et all Il997t iDeharveng et all l200ll) . it is pos- 
sible to evaluate the observed lower limits of -Fuv/-Flc 
for comparison with higher z objects (the bo ttom part of 
Table 4; Table 1 in IDeharveng et all l200l|) . Except for 
Mrk 54, these limits, corrected for the foreground absorp- 
tion by the Galactic HI and H2 gas, are not better than 
those obtained at high- 2: (column 3). The low luminos- 
ity of the sample galaxies is one of the causes (column 
7). Nevertheless, moderate upper limits are obtained for 
the relative escape fraction, < 25-40% (3-cr confidence 
level; column 5) , because no correction for the IGM opac- 
ity is necessary. For Mrk 54, whose luminosity is similar to 
th ose of our sample (colum n 7), the observed lower limit 
of IDeharveng et aD |2001) translates into an upper limit 
on the relative escape fraction of 3% comparable to those 
for z ~ 1 galaxies of M03 (column 5). 

5.2. Absolute escape fraction 

As discussed above and shown by equation (2), estimat- 
ing the absolute escape fraction from the relative escape 
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fraction requires an evaluation of the dust attenuation 
withi n each galaxy which is often difficult and uncertain 
(e.g iBuat et all l2002t iPettini et aflll998llMeurer et all 
Il999l) . Although the infrared to U V flux ratio is a 
good estimator of t he UV attenuation ijBuat et allll999t 
iGordon et afll2000l) , the infrared fluxes for high-z galaxies 
are not available yet. Here, we adopt a calibration between 
the noni onizing UV sl o pe an d the UV attenuation pro- 
posed bv lMeurer et alJ l|l999t) for simplicity. However, we 
should keep in mind that the cali bration depends on the 
type of galaxies, starburst or not l)Bel]Ll2002tlKong et all 
120041) . 

The UV slope of the galaxies listed in Table 4 have 
been searched in the literature, or, if not available, esti- 
mated from broad-band measurements with the assump- 
tion of a power-low spectrum (f\ tx A' 5 ) 6 (references in 
column 9 of Table 4). This was not possible for some 
galaxies of M03 because their colors in the rest-frame 
A < 3000 A are not available. For the composite spec- 
trum of SOI, we estimated the UV attenuation from the 
repo rted mean E(B — V ) via the Calzetti's attenuation 
law ((Calzetti et alll200C|) . The UV slope of the composite 
spectrum shown in Fig.l of SOI is consistent with the slope 
corresponding to the estimated attenuation. We estimated 
UV slopes of our two galaxies from their broad-band colors 
although we have UV spectra of our two galaxies because 
the data quality is not good. 7 

Since the calibration gives the attenuation at the rest- 
frame 1600 A, we convert it into the attenuation at the 
app ropriate UV waveleng th by the Calzetti's attenuation 
law ijCalzetti et allEoOOli if the UV wavelength is different 
from 1600 A (column 8). We note here that the uncertainty 
resulting from those on the UV slope and colors is very 
large; for example, the uncertainty of about 0.05 mag in 
^606 — -^814 for our two galaxies translates into AAttv ~ 
2A0 ~ 6A(V r 606 - Jsu) ~ 0.3 mag (see iMeurer et all 
119991) . 

For galaxies at z ~ 3, we find that (1) the absolute 
escape fraction from the detection by S01 is 60 ± 13% 
(l-cr observational uncertainty), (2) the absolute escape 
fraction of the two brightest LBGs observed by G02 is 
less than 5% (3-cr), (3) the absolute escape fraction of L* 
LBGs observed by us is less than 20-40% (3-cr). 

For z ~ 1 galaxies, we find very small upper limits on 
the absolute escape fraction, typically, less than a few per- 
cent (3-<r). Since small upper limits were obtained even for 
the relative escape fraction, the conclusion of very small 
absolute escape fractions for the observed galaxies seems 
robust against the uncertainty of estimating the dust at- 
tenuation. 



In the estimation, we neglected the effect of the IGM HI 
absorption on the broad-band photometry because it is small 
enough, although this results in an overestimation of /3. 

7 We can measure the UV slopes in the spectra. The obtained 
slopes for the two objects are similar, (5 = —1.9 ± 0.3, which 
translates into Auv = 0.7 ± 0.6 mag. 



For nearby galaxies, we find small upper limits on 
the absolute escape fraction, less than 10% (3-cr). These 
values are a factor of 0.1-1 smaller than those esti- 
mated from t he comparison between the LC fluxes and 
the Ha fluxes llLeitherer et allll993lHurwitz et allll997t 
iDeharveng et allfeOOllT T cause for this discrepancy may 
be the LC extinction by dust in H 11 regions; LC photons 
are absorbed b y dust before they ionize neutral hydro- 
gen a toms fe.g.. llnoue. Hirashita. fc Kamavail200lUlnoueL 
200 lj). This effect leads us to underestimate the intrinsic 
LC flux from the flux of the recombination line and to 
overestimate the escape fraction. However, uncertainties 
of UV and Ha attenuations are also likely to play a role. 

6. Conclusions 

We made the first attempt to constrain the amount 
of Lyman continuum (LC) escaping from galaxies, us- 
ing narrow-band photometric observations with the 
VLT/FORS, and then, reached the following conclusions: 

1. Because of an unexpected systematic effect on pho- 
tometric redshifts, only two objects with a redshift 
appropriate for measuring the LC radiation are left 
in our sample. None of these two galaxies, HDFS 85 
(z sp = 3.170) and HDFS 1825 (z sp = 3.275), which are 
~ L* Lyman break galaxies (LBGs), are detected at 
a 3-cr level. The resulting 3-cr lower limits of the ob- 
served UV-to-LC flux density (per Hz) ratio are 15.6 
and 10.2 for HDFS 85 and 1825, respectively. These 
limits are compatible with the detection in the com - 
posite spectrum of 29 LBGs bv ISteidel et~afl l|200lh . 
iGiallongo et"afl (2002) obtained slightly better limits 
than ours because of the high luminosity of their two 
galaxies. 

2. After a comparison with population synthesis models 
and a correction for average IGM opacity, the observed 
lower limits translate into relative escape fractions less 
than 0.7 and 2.2 (3-cr confidence level) for HDFS 85 
and 1825, respectively. 

3. In addition to the two objects observed by 
IGiallongo e~ al. ( 2002), HDFS 85 makes the third case 
of a relative escape fra ction smaller than that reported 
bv lSteidel et all |200l|) . Fluctuations of the IGM opac- 
ity from line of sight to line of sight and randomness 
of LC escape from object to object may play a role in 
this discrepancy. 

4. Very small values of the absolute escape fraction are 
estimated from the relative escape fraction, namely less 
than 10%, except for a few cases. The dust attenuation 
which is difficult to evaluate from the UV data only is 
a source of uncertainty. 

5. When the high l umin osity of the galaxies observed by 
IGiall ongo et al. (2002) is accounted for, our observa- 
tions show that narrow-band photometry can reach 
stronger limit than spectroscopy in terms of the rel- 
ative escape fraction. This advantage is not obtained 
with broad-band imaging at high-z because of the IGM 
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opacity. Stacking a significant number of deep narrow- 
band images of drop-out galaxies has, therefore, the 
potential to confirm or not the high re lative escape 
fraction reported bv ISteidel et alJ l)200l|) . In addition 
to increasing sensitivity, such a method would average 
the IGM opacity and the randomness of the LC escape. 
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Appendix A: A mean IGM opacity model 

In this appendix, we first show that the number density distribution function of the Lyman limit systems (LLSs) 
reasonably agrees with an extension of the distribution function of the Lyman a forest, in contrast to previous works. 
That is, the intergalactic clouds may be described by a continuous distribution from low- to high-density through any 
redshift. Second, we present an analytic approximation of a mean opacity of the intergalactic medium (IGM) based 
on the distribution function. 

A mean IGM opacity at an observ ed wavelength A Q bs against a source with z = zg can be expressed as (e.g., 



A mean lUlvi opacity at an observe* 
iMrtller fc Ja,kob^ir^lMadauiri99i 



rzs rNm,*j> f) 2 \f 

(t{™(z s ))= dz dNm ^J^[ 1 - exp{ - (JH (X ohs /l + z)N m }], (A.l) 



o «/JV H i,iow dzdN m 

where Njjx is the HI column density of a cloud in the IGM, A^htjow and A^Hi.up are the lower and upper limit of the 
cloud distribution, and cth(A) is the hydrogen cross section at the wavelength A. Here, we assume a differential number 
distribution function of the IGM clouds expressed as 

f(Nm)g(z) , (A.2) 



dzdN m 



where /(Ahi) oc N^f with the normalization as J^™'^ v f(Nm)dN-Hi = 1 and g(z) = A{\ + z) 7 . We note here that 
the normalization A means the total number of the IGM cloud per unit redshift interval and is a function of A^hi.Iow 
and A^HLup- 

Observationally, the power-law index (3 may be regarded as a constant equal to about 1.5 for a wide range of 
colum n densities from log Am /cm- 2 = 12.5 to 21 5, althoug h it se ems to depend on re dshift and column density in 
detail ((Sargent et all Il989t IPetitiean et all Il993t [Kim et all l2002t [Peroux et all l2008h . For simplicity, we assume a 
single power-law index (3 — 1.5 for all redshifts a nd column den sities. 

For the cloud number evolution with redshift. iMadaiJ jl995h adopted two different evolutionary in dices (7) for the 
Lyma n a f orest {Nm ^ 1 17 cm " 2 ) and for the LLSs (A^hi iS 10 17 cm -2 ), based on observations bv iMurdoch et all 
l)l986|) and ISargent et aD l[l989). However, the same index 7 for these populati ons is compatible w ith recent ob- 
servations as shown in Fig. A.l. The LLS number evolution is reproduced from IPeroux et aD J2003) in this figure 
as the filled square point s with error bars . The solid line is the number evolution of the Lyman a forest reported 
by iKim et aD (2002) and Wcvmann et al multiplied by the reduction factor due to the lowest column den- 

sity co nsideredT lKim et alJ l)2002|) deal with Lyman a forest with logAfHi/cm~ 2 = 13.644-17, whereas IPeroux et alJ 
l)2003|) deal with the clouds with a column density larger than 1.6 x 10 17 cm~ 2 , so that the multiplicative factor is 
(4.4 x 10 13 /1.6 x lO 17 )^ 1 « 0.017 when (3 = 1.5. There is a systematic disagreement along the ve rtical axis for 
high-z points between the solid line and the data, but the evolutionary slopes are consistent as noted bv lPeroux et alJ 
(2003). The high-z data are repro duced by the dashed line whose redu ction fa c tor is 0.01. On the other hand, the 
LLS number evolution reported bv lSargent et al.l (1989) and adopted by iMadaul (1995) is shown as the dotted line in 
Fig. Al and cannot r e produ ce the data points of IPeroux et alJ (|2003|) . Although the number evolution suggested by 
IStengler-Larrea et al.l l)l995[) shows a better fit (dash-dotted line), in this paper, we assume the case of the solid line. 
This means that we assume the same redshift evolution of the cloud number density over the entire column density 
range with (3 = 1.5. 

Under a cloud distribution function with a single power-law index (3 (1 < (3 < 2) and the same evolutionary index 
7 for all range of the column density, we take a set of limiting column densities as o'hA'hIjIow *C 1 and (ThA^hi^ip 3> 1, 
for example, A^hliow = 10 12 cm~ 2 and A^Hi, U p = 10 22 cm~ 2 . In this case, equation (A.l) can be approximated as 

MSf (*s)> « / r(2 - (3)AN£^ W (1 + z)V H (A obs /l + zf-'dz , (A.3) 
Jo 

where r(2 — f3) is the usual Gamma function l|Zuo & PhinnevL Il993h . The factor, ^lA^^j^ can be estimated from 
the observed number of clouds for a limited range of the column density with the power-law distribution. If the 
number of clouds with a column density between N\ and N u (Ahi,1ow < A) < Au < A^ni.up) is denoted as A', we 
have ANhtw ~ A'Nf ' 1 when A^hlup/A^hijow > 1, Ar u /AT H i,iow > 1, and (3 > I. According to lKim et all l|2002l) and 



IWevmann et all l(l998j) . we adopt (^',7) = (6,2.5) for z > 1.1 and (34,0.2) for z < 1.1 against log Am/cm- 2 = 13.64- 
17 (i.e. Ni = 4.4 x 10 13 cm" 2 ). 

We will present an analytic approximation of equation (A.3), although we can integrate it numerically with a 
detailed function of the hydrogen cross section, cth- We adopt an approximated form of the cross section for the 
analytical formula. The point is that the cross section of the i-th line, <7i(A), can be neglected for a wavelength A out 
of a small range of |A — A,| ^ b/c, where Aj is the central wavelength of the line, b is the Doppler parameter, and c is 
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Fig. A.l. Number density evolution of the Lyman limit systems (LLSs; l ogA^HT/cm" 2 = 17.2 -20) in a unit redshift 
interval along a line of sight. Data points with error-bars are taken from IPeroux et all l|2003|) . Th e solid line is the 
nu mber evolution expected from that of the Lyman a forest (log iVm/cm -2 = 13.64-17) reported bv lKim et all ll2002h 
and lWevmann et alJ 1 199S|) with a single power-law distribution of the column density; the number of LLSs is that of 
Lyman a forest multiplied by (4.4 x 10 13 /1.6 x 10 17 )^ -1 w 0.017 when (3 — 1.5. The dashed line is th e same as the 
solid l ine but the multiple factor is 0.01. The dotte d line is the number evolutio n of LLSs reported bv lSargent et alJ 
l)l989|) . and the dash-dotted line is that reported bv lStengler-Larrea et al.l l|l995|) . 



the speed of light. Besides, the cross section for LC, ctlc(A), is zero when A > Al = 912A. Thus, we can approximate 
°h(A)' 3 ~ 1 (= [olc(A) + Yli <J i(A)]' 3_1 ) t° <?lc(A)^ _1 + J2i cr i(^)' 3 ~ 1 ! because only one term has meaningful value for a 
specific wavelength A. 

Adopting this approximation, we can reduce equation (A. 3) to 

where 



and 



, \ 303-1) 

T(2 - /3)(a L iV H i,iow) /3 - 1 (^) ' J^ m A(l + zy-W+Hz (for A obs < A L (1 + z s ) ) ^ (A 5) 

(otherwise) 

/3)(a , l 7V H i,iow) /3 -U(^) (^xf) 7+1 (forA l <A obs <A l (l + z s )) , (A6) 

(otherwise) 

For the Lyman continuum, <th(A) = ol(A/Al) 3 with ox = 6.30 x 10 -18 cm 2 ijOsterbrockl 179891 being the cross section 
at the Lyman limit. The integral in the last term of equation (A. 5) can be done analytically if z dependences of 7 
and A are simple. The lower edge of the integration is z m - ln — max[0, (A OD s/Al) — 1]- For Lyman series lines, whose 
central wavelength and cross section are denoted as Xi and <7o,i, we have integrated the integral in equation (A. 3) over 
a narrow range characterized by the Doppler parameter, 6, around the line center, assuming a rectangle line profile. 
Actually, the integration has been made over the range of the Doppler width multiplied by a numerical factor 8 which 
is order of unity. From the comparison between a numerical integration of equation (A. 3) with a detailed function of 
hydrogen cross section 8 and the analytic approximation of equation (A. 4-6), we find that the appr oximation become s 
excellent if S = 2.6 for b — 30 km s^ 1 , which is a typical Doppler parameter for Lyman a forest l|Kim et all l2002|) . 

8 We approximately take into account the Voigt profile by the following equation; (f>(x) ~ exp(— x 2 ) + a/^Tvx 2 , where a is the 
damping constant, and x is the relative displacement in the frequency space, i.e. x = (u — i>i) /Avd, where v, Vj, and Az^p are 
the frequency, the central frequency of the line, and the Doppler width in the frequency space, respectively jFerlandl Il996'). 
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Fig. A. 2. Examples of the IGM optical depth calculated by an analy tic approximation presented as equations (A. 4-6). 
The dotted curve is the model "MC-Kim" bv lBershadv et al 1 lll999h . 



All c alculations include up to 40-th line of Lyman series whose central cross sections are calculated by the data taken 
from lWiese et alJ l|l96fih . 

Fig. A. 2 shows the mean I GM opacit i es for several source redshifts calculated by equations (A. 4-6) with the 
distribution function based on iKim et al.l ()2002t) and (5,b) — (2.6, 30kms _1 ). Since Kim's observations are made 
against Lyman a forest below z ~ 4, we have sh o wn ca ses up to z$ — 4. As a comparison with a previous opacity 
model, we also plot the result bv iBershadv et all lll999l) fthe case of "M C-Kim") as the dotted cu rve in the figure. 
Since they adopted the LLS number evolution of ISargent et al.1 l|l989(l as done bv iMadaul l|l995l) . there are some 
discrepancies between ours and theirs; our larger (smaller) optical depth at around the obser ved 3000 (150 ) A th an 
IBershadv et alJ l|l999f) is caused by our larger (smaller) number of LLSs at z ^ 2 (^ 2) than lSargent et all ljl989|) as 
shown in Fig. A.l. However, the quantitative agreement is reasonably good. 

As another test, we compare our optical depth with that measured bv ISteidel et alJ l 200lh. From a comp arison 
between two composite spectra of galaxies and QSOs at (z) ~ 3.4 and 3.3, respectively. ISteidel et all HflOl) have 
found a mean IGM optical depth of 1.35 for the rest-frame 900 A, whereas our model gives 1.36 for the wavelength 
and z = 3.4 source. We have found an excellent agreement with each other. 



